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Abstract 

The nineteen-vertex models of Zalomodchikov-Fateev, Izergin-Korepin and 
the supersymmetric osp(l|2) with periodic boundary conditions are studied. 
We find the spectrum of these quantum spin chain using the Coordinate Bethe 
Ansatz. The approach is a suitable parametrization of their wavefunctions. We 
also applied the Algebraic Bethe Ansatz in order to obtain the eigenvalues and 
eigenvectors of the corresponding transfer matrices. 
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1 Introduction 



One-dimensional quantum spin chain Hamiltonians and classical statistical systems 
in two spatial dimensions on a lattice (vertex models), share a common mathemat- 
ical structure responsible by our understanding of these integrable models [|l], 0. If 
the Boltzmann weights underlying the vertex models are obtained from solutions of 
the Yang-Baxter (YB) equation the commutativity of the associated transfer matrices 
immediately follow, leading to their integrability. 

The Bethe Ansatz (BA) is the powerful method in the analysis of integrable quan- 
tum models. There are several versions: Coordinate BA [^, Algebraic BA Q, Analyt- 
ical BA [^], etc. developed for diagonalization of the corresponding Hamiltonian. 

The simplest version is the Coordinate BA. In this framework one can obtain the 
eigenfunctions and the spectrum of the Hamiltonian from its eigenvalue problem. It is 
really simple and clear for the two-state models like the six-vertex models but becomes 
awkward for models with a higher number of states. 

The Algebraic BA, also proverbial as Quantum Inverse Scattering method, is an 
elegant and important generalization of the Coordinate BA. It is based on the idea of 
constructing eigenfunctions of the Hamiltonian via creation and annihilation operators 
acting on a reference state. Here one uses the miraculous fact that the YB equation 
can be recast in the form of commutation relations for the matrix elements of the mon- 
odromy matrix which play the role of creation and annihilation operators. From this 
monodromy matrix we get the transfer matrix which, by construction, commutes with 
the Hamiltonian. Thus, constructing eigenfunctions of the transfer matrix determines 
the eigenfunctions of the Hamiltonian. 

Imposing appropriate boundary conditions the BA method leads to a system of 
equations, the BA equations, which are useful in the thermodynamic limit. The energy 
of the ground state and its excitations, velocity of sound, etc., may be calculated in 
this limit. Moreover, in recent years we witnessed another very fruitful connection 
between the BA method and conformal field theory. Using the Algebraic BA, Ko- 
repin found various representations of correlators in integrable models and more 
recently Babujian and Flume developed a method from the Algebraic BA which 
reveals a link to the Gaudin model and render in the quasiclassical limit solutions of 
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the Knizhnik-Zamolodchikov equations for the SU{2) Wess-Zumino-Novikov-Witten 
conformal theory. 

Integrable quantum systems containing Fermi fields have been attracting increasing 
interest due to their potential applications in condensed matter physics. The prototyp- 
ical examples of such systems are the supersymmetric generalizations of the Hubbard 
and t-J models [§. They lead to a generalization of the YB equation associated with 
the introduction of the a Z2 grading |^ which leads to appearance of additional signs 
in the YB equation. 

In this paper we consider the Coordinate and Algebraic versions of the BA for the 
trigonometric three-state vertex models of 19-vertices with periodic boundary condi- 
tions. These models are well-known in the literature: the Zamolodchikov-Fateev (ZF) 
model or A\ model [0, the Izergin-Korepin (IK) model or model [ill] 
supersymmetric osp{l\2) model [p!2| . 



While the BA solution of the periodic ZF model was derived by a fusion procedure 
in ||T^ and [|r3|, a generalization of the Algebraic BA was developed by Tarasov 
16| to solve the IK model. Moreover, the IK model was solved via Coordinate BA by 



Batchelor at al in [17]. 

In the context of the Algebraic BA , the version presented here is based on the 
Tarosov approach but now we include the ZF model and we also extend it to the 
graded version of the quantum inverse scattering method in order to consider the 
osp{l\2) model. 

In the context of the Coordinate BA, we propose here a new parametrization of 
wavefunctions. This result is of fundamental importance since it allows us to treat 
these 19-vertex models in the same way and the Coordinate BA for these three-states 
models becomes simple enough as for two-state models. 

The main goal in this paper is to reveal the common structure of these 19-vertex 
models which permits us to apply the BA method, unifying old and new results. 

The paper is organized as follows: In Section 2 we present the models. In Section 
3 the spectra of the corresponding Hamiltonians are derived using the Coordinate 
BA and in Section 4 the Algebraic BA is also used to diagonalize the corresponding 
transfer matrices. We justify this twofold presentation remarking that the BA method 
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is apparently version dependent. It means that when one solves a model using a 
particular BA version is not always clear how to extend the solution for all possible 
versions. For example, the biquadratic model was solved by Coordinate BA in [0, |T^ 
and its Algebraic BA version is still unknown. Finally, the conclusions are reserved for 
Section 5. 

2 Description of the models 

Let us start with the graded formulation and then recover the non graded from it. 

Consider V = Vq Q) Vi a. Z2-graded vector space where and 1 denote the even 
and odd parts respectively. The multiplication rules in the graded tensor product 

s 

space V ^ V differ from the ordinary ones by the appearance of additional signs. The 

s s 

components of a linear operator A ^ B & V ^ V result in matrix elements of the form 

{A ® 5)5 = (_)p(/3){pH+i>(7)) a^^B^s (2.1) 

s s 

The action of the graded permutation operator V on the vector \a) ® \P) & V ® V is 
defined by 

V \a) ® 1/3) = (-)^'W^'(^) 1/3) ® \a) =^ {V)tp = {^Y^^^p^^^ 6^s (2.2) 

where p{q) = 1 (0) if \a) is an odd (even) element. 

The central object in the theory of integrable models is the 7?.-matrix TZ{X), where 
A is the spectral parameter. It acts on the tensor product V'^ (S> V"^ for a given vector 
space V and it is solution of the Yang-Baxter (YB) equation 

7^12(A)7^13(A + /i)7^23(;u) = 7^23(/i)7^l3(A + /i)7^l2(A) (2.3) 

in (g)V^ (g) \/^ where 11x2=11® X, 7^23 = X ® 7^, etc. 

In the graded case, 7?.i3 however, does not act trivially on the second space due to 
signs generated by commuting odd operators. In this case, the graded YB equation in 
components reads 

^^''(A)7^£^(A + ^)<,{:(/.)(-)^(^')(^(^")+^(^")) 

= 7^yf (/i)7^^^■:(A + /i)7^^{,(A)(-)^('=')(^(^■")+^('=")) (2.4) 
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Besides TZ we have to consider matrices R = VIZ which satisfy 

i?12(A)/?23(A + li)RM = i?23(/i)i?12(A + /i)i?23(A) (2.5) 

Because only R12 and -R23 are involved, Eg . (12 .51) written in components looks the same 
as in the non graded case. Moreover, the matrices TZng = PR satisfy the ordinary YB 
equation (|2.3| ) where P is the non graded permutation operator. 

2.1 The i?-matrices 



We will consider 19-vertex models for which their R matrices have a common form 
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(2.6) 



Here we have assumed that the grading of threefold space is p(l) = 1, p{2) = e and 
p{3) = 1, where e = ±. The matrix elements Xi and r/i for each model will be listed 
below . 

• ZF R-matrix 

The simplest 19-vertex model is the ZF model or A\ model . The solution of the 
YB equation was found in ||10[. It can also be constructed from the six- vertex model 
using the fusion procedure. The corresponding i?-matrix has the form ( |2.6| ) with e = 1 
and 



Xi(A) = sinh(A + 77) sinh(A + 2?7), X2(A) = sinh A sinh(A + r^) 

X2,{X) = sinh A sinh(A — r]), X5(A) = y5(A) = sinh(A + 77) sinh 27] 

xq{X) = ye{X) = sinhAsinh2?7, X7(A) = yji^X) = sinh 17 sinh 27] 

X4(A) = X2(A)+X7(A) 



(2.7) 
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IK R-matrix 



The solution of this YB equation was found in [|ri[] . It can not be constructed from 
the six-vertex model using the fusion procedure. The i?-matrix has the form ( |2.6|) with 
e = 1 and 



sinh2r7(e^"3'? + e3^) 



sinh(A — 5i]) + sinh?], a;2(A) = sinh(A — 3r]) + sinhSr] 
sinh(A — rj) + sinhr] 

sinh(A — 31]) + sinh Si] — sinh brj + sinh rj 
- sinh 2r/(e-^+3'' + e'^'^), y^iX) = 
e2'^sinh2r7(l -e"^), ye^X) = e'^^ sinh 2r/(l - e^) 
— 2e~'^~''^'' sinh r] sinh 2r] — sinh Arj 
2e^~^^ sinh rj sinh 2?] — sinh Arj 
• Osp{l\2) R-matrix 

The trigonometric solution of the graded YB equation for the fundamental repre- 



xi(A) 
xsiX) 
Xi{X) 

X5(A) 

X(i{X) 
xjiX) 
yriX) 



(2i 



sentation of osp{l\2) algebra was found by Bazhanov and Shadrikov in [|T^]. It has the 
form (|2.6| ) with e = —1 and 

sinh(A + 2ri) sinh(A + Srj), X2{X) = sinhAsinh(A + Srj) 
sinh A sinh(A + rj) 

sinh A sinh(A + Si]) — sinh 2// sinh 3rj 

e"'^/^ sinh2?7sinh(A + Sr]), ?/5(A) = e''*/^ sinh 2?] sinh(A + Sr^) 
_Q-x/3~2r, gj^j^ sinh A, ye (A) = e^/^+^'' sinh 2?? sinh A 
e^/^ sinh 2?7 (sinh(A + 3r]) + e"'' sinh A) 
e"^/^ sinh 2?] (sinh(A + 3r]) + sinh A) 



Xi{X) 

X3(A) 
X4(A) 
X5(A) 
Xq{X) 

x-ji^X) 
y7(A) 



(2.9) 



The rational limit of ( |2.9| ) is well-known in the literature PD| and can be written 
in the form : 

Xt] 



R{X,r]) =r]I+XV 



A + f 



(2.10) 



where T is the identity operator, V is the graded permutation operator ( |2.2| ) and U 
is the rank-one projector = U. The algebraic solution of ( |2.1CI| ) was obtained by 
Martins |21|, as a limit of the algebraic solution of the IK model. 
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2.2 The Hamiltonians 



In order to derive the Hamiltonian, it is convenient to expand the i?-matrix around 
the regular point A = 0. For the 19-vertex models the corresponding solutions with 
the standard normalization can be read directly from ( |2.6| ). They have the form 

(2.11) 



R{X, ri)r^I+ \{a-^n + pi) + o(A2). 



with a and P being scalar functions. 

The Hamiltonian is a local sum given by 



k=l 



(2.12) 



where H^^k+i is the H in (p.ll|) acting on the quantum spaces at sites k and k + 1. Using 
a spin language, this is a spin 1 Hamiltonian. In the basis where 5*^ is diagonal with 
eigenvectors |+, k) , |0, k) ,\—, k) and eigenvalues 1, 0, —1, respectively, the Hamiltonian 
densities acting on two neighboring sites are then given by 
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(2.13) 



k,k+l 

where the matrix elements for each model are: 
• ZF Hamiltonian 

For the ZF model the corresponding quantum spin chain is the spin 1 XXZ model. 
The two site Hamiltonian is derived from ( p.ll| ) and has the form ( p.l3| ) with 



e 

Zl 
Z5 
Z7 



1, a = sinh 2?7, /5 = 
0, 2:3 = —1, Z4 = —2 cosh2ri, 
z^= — cosh 2?7, zq =Zq= 2 cosh r] 
Zj= —1 — 2 cosh2?7 



(2.14) 



6 



• IK Hamiltonian 



In the IK model the two-site Hamiltonian for the corresponding quantum chain has 
the form (|2.13| ) with 



£ = 1, a = — 2sinh2r7, /? = 

cosh f] sinh At] sinh rj 

= 0, = — — — , Zi = —2 — , 

cosh Sr] cosh drj 

2„ sinh 2rj _ _^ sinh 2?7 

^6 = ^ — \ — ^6= -e ^ — - — - 
cosh 37] cosh 3?7 

^ ^-4, coshr/ ^ ^4. coshr^ 

cosh 3?7 ' cosh 3?7 

• Osp(l|2) Hamiltonian 

The two site quantum Hamiltonian associated with the osp(l|2) model has the form 
with 



e = —1, a = sinh2?7, /? = — coth2?7 

sinh?7 

zi = cosh277, zs = — , ^4 = 1 + coth3r7smh277, 

sinh 3r] 

sinh 2?7 

Z5 = ^ — 5 ^5= —^5 

n„ sinh 2?7 _ sinh 2r7 

^6 = -e ^ — - — -, ^6= e — - — - 
sinh 3?7 sinh 3r] 

^ sinh2r/ ^ sinh 2?/ ^ ^ _ sinh2r/ ^ ^^ sinh2r^ ^ 
3 sinh3r7' 3 sinh3r7 

Having now built a common ground for these models, we may proceed to find their 
spectra. We begin with the Coordinate BA because of its simplicity. 



3 The Coordinate Bethe Ansatz 

In this section results are presented for a periodic quantum spin chain of atoms each 
with spin 1 described by the Hamiltonian ( |2.12| ). At each site, the spin variable may 



be +1, 0, —1, so that the Hilbert space of the spin chain is H^^^ = where V = 
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with basis {|+) , |0) , |— )}• The dimension of the Hilbert space is diniTi^^^ = 3^. On 
7-^(^) we consider the Hamiltonians presented in the previous section. 

From ( |2.12| ) one can see that H commutes with the operator which shifts the states 
of the chain by one unity. It means translational invariance of H. Moreover, the 
Hamiltonian ( |2.12|) preserves the third component of the spin 

N 

[^,^t] = 0, ^T = E^fc (3-1) 

k=l 

This allows us to divide the Hilbert space of states into different sectors, each labelled 
by the eigenvalue of the operator number r = N — . We shall denote by 7i^^^ the 
subspace of 7i^^^ with r = n. One can easily see that dimTi*^^^ = Y,^=o dimH^J^^ with 

where [|] means the integer part of |. 

3.1 Sector r=0 

The sector Ti.Q^'' contains only one state, the reference state, with all spin value equal 
to +1, \I^o = 11^1+,^), satisfying H'^o = -Eq^o, with Eq = Nzi. All other energies 
will be measured relative to this state. It means that we will seek eigenstates of H 
satisfying {H — Nzi)'ifr = e^^l/r , in every sector r. 

3.2 Sector r=l 

In 'Hi^\ the subspace of states with all spin value equal to +1 except one with value 

+ + +0 + + ■ — h) which span a basis of Hi^\ The 

k / 



0. There are N states \k[0]) 



ansatz for the eigenstate is thus of the form 

N 

^, = j:Aik)m) (3.3) 

k=l 

The unknown wavefunction A{k) determines the probability that the spin variable has 
the value at the site k. 

From the complete invariance translational due to the periodic boundary conditions, 
it follows that A{k) is just the wavefunction for a plane wave 

Aik) = e (3.4) 
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where ^ = e^^ , 6 being some particular momentum fixed by the boundary condition 
A{N + k) = A{k). 

When H acts on |A;[0]) , it sees the reference configuration, except in the vicinity 
of k, and using (|2.13|) we obtain the eigenvalue equations 

(ei + 2zi - Z5- Z5)A{k) = A{k - 1) + A{k + 1) (3.5) 

The plane wave parametrization ( p. 41) solves (|3.5|) provided 

ei = -2zi + Z5+ ^5 +e + r' (3.6) 

Thus \E'i is the eigenstate of H in the sector r = 1 with eigenvalue Ei = {N — 2zi) + 
Z5+ ^5 +2 COS 6*, where 9 = 2'kI/N ,1 = 0,1, ...,N - 1. 

3.3 Sector r=2 



In the Hilbert space Ti^^^ we have states of the type ]) 

and N{N -l)/2 states of the type |A;i[0], A;2[0]) 
these eigenstates in the form 



+ +-+ + ■■■ + 

k 

+ + + + + H h y We seek 

ki k2 I 



N 

^2= E Aik,,k2)\k,[0],k2m + J2B{k)\k[-]). (3.7) 

ki<k2 k=l 

The periodicity condition reads now 

A{k2,N + ki) = A{ki,k2) and B{N + k)=B{k) (3.8) 

Following Bethe the wavefunction A{ki, ^2) can be parametrized using the superpo- 
sition of plane waves (|3.4|) including the scattering of two pseudopartides with momenta 
e, and 62, (e, = e^^^j = l,2): 

Aik,, k2) = A,2ei'i2' + ^le^e'^ (3.9) 
which satisfy the periodic boundary condition ( p.8|) provided 

A2 = ^2ief , ^21=^2^2"" (3.10) 

and the parametrization of B{k) is still undetermined at this stage. 
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Before we try to parametrize B{k) let us consider the Schrodinger equation {E2 — 
Nzi)'$2 = £2^1/2 • From the exphcit form of H acting on two sites ( |2.13| ) we derive the 
following set of eigenvalue equations 

(e2 + 4zi - 22:5 - 2 55)A(/ci, ^2) = A{ki-l,k2) + A{ki + l,k2) 

+A{k^,k2~l) + A{k,,k2 + 1) (3.11) 
{e2 + 3zi- z:^- Z5 -ez4)A{k,k + l) = A{k - l,k + 1) + A{k,k + 2) 

+e Zq B{k + 1) +€ZGB{k) (3.12) 
{e2 + 2zi- Z7- Zr)B{k) = z^Bik - 1) + z^Bik + 1) 

+ Zq A{k -l,k) + zeA^k, k + 1) (3.13) 

The parametrization ( p.9|) solves the equations ( p.ll|) provided 

62 = -Azi + 22:5 + 2 z, +^1 + + 6 + ^2"' (3-14) 

It immediately follows that the eigenvalues of H are a sum of single pseudoparticle 
energies. 

The parametrization of B{k) can now be determined in the following way: Sub- 
tracting Eq. (|3.12|) from Eq. (|3.11|) for ki = k, k2 = k + 1, we get a meeting equation 



e z^ B{k + l)+6Z6B{k) = A{k, k)+A{k+l, k + l)-{zi+ezi-Z5- Z5)A{k, k + l) (3.15) 

Now we extend the parametrization ( p.9| ) to ki = k2 in order to get a parametrization 
for the wavef unction B{k) : 

B{k) = 5(^16)' (3.16) 
which solves the meeting equation ( |3.15| ) provided 

1+^16- A16 , ^ 1+^16-Ai^. 

Z6+ 4l42 Zq+ Zq ^1^2 

Ai = zi+ ezi - Z5- Z5 (3.17) 

These relations tell us that the pseudoparticle of the type \k[—]) behaves as the two 
pseudoparticles |fci[0]) and |A;2[0]) at the same site k and its parametrization follows as 
the plane waves of particles \ki[0]) multiplied by the weight function B = B{C,i,^2)- 
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Now substituting ( pl9| ) , ( p.l4| ) and ( |3.16| ) into the equation ( p.l3| ) we find the phase 
shift of two pseudoparticles. 

(1 + - (1 + 66)(A26 + A36) + + As^l 

A,2 (l+66)'-(l+66)(A2e2 + A3ei)+A4eie2 + A5e? ^' ' 

where 

A. ^ 1 

A3 = h —(2:32:4 - Zq Zq) + {Zi - Z^- Z^) 

z^ 

A4 = —{ezi + z^+z■J)^ {zi- zr^ -z^) -2 

Z3 

A5 = —{ez^ + zi-z^-z^) (3.19) 
Z3 



Combining this resuh with the periodic relations ( p.lO|) and using (|2.14| - ^TT6D we arrive 
to the BA equations in 7^2^^ for each modeh 



.N / 1 + 66 + 6 + 6 - (A + 2)^2 



~ U + a2 + 6 + 6-(A + 2)6j ^^-^^^ 



for the ZF model, 



^ ^ _ / i + ei6-A6 \ A + 66 -6 -6 + (A- 2)^1 \ 

U + 66 - AG Ai + 66 - 6 - 6 + (A - 2)67 ^ ' ^ 

for the IK model and for the osp{l\2) model we get 

_ / i + ei6-A6 \ A + 66 + 6 + 6- (A + 2)^1 \ 
' lvi + 66-AeiAi + 66 + 6 + 6-(A + 2)6y ^' ^ 

where 

A = 2cosh2r7, and {^^^2)^ = 1. (3.23) 

3.4 Sector r=3 

Now the Hilbert space is 7^3^^ where there are A^(A^ — 1)(A^ — 2)/6 states of the type 
\ki[0], k2[0],k3[0]), N{N - l)/2 states of the type |A;i[-], A;2[0]) and N{N - l)/2 states 
of the type |fci[0], A;2[— ]) . We seek eigenfunctions of the form 



^3 = E A{k,,k2,k,)\k,[0],k2[0],ksl 
ki<k2<k3 

+ E {Bliku k2) k^m + fi2(A:i, k^) |A;i[0], k^i-])] (3.24) 



k<k 
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Periodic boundary conditions read now 

A{k2, h, N + h) = A{h, k2, h), B^ih, N + h) = B,{h, k2), (3.25) 

Again, the wavefunctions A{ki,k2,k3) can be parametrized by the superposition of 
plane waves 

A{k,, k2, k,) = Au3ei'^2'& + AlS2ei'^2'& + A2l3e^^2'& 

+A23iei'e2'^t + As^ei'^t^t + ^123^1^e2^e3^ (3.26) 

which satisfy the periodic boundary condition provided 

^231 _ ^321 _ ^213 _ ^312 _ ^123 _ ^132 _ fn nj\ 

4 ~4 -^3) 1 -1 - ^2 > -1 - -1 - ^1 [6.Z{) 
^123 ^213 ^132 ^231 ^312 ^321 

These relations tell us that the interchange of two pseudoparticles is independent of 
the position of the third pseudoparticle. Using 5'-matrix language, this locality of the 
interactions is equivalent to the factorization property of the S'-matrix, according to 
which the scattering amphtude of three particles factorizes into a product of three 
two-point S-matrices. 

Action of H on these eigenstates gives the following set of coupled equations for 
A{ki, k2, ks) and Bi{ki, ^2), i = 1, 2: 

(eg + 6zi - 3^5 - 3 Z5)A{ki, fca, ks) 
= A{ki-l,k2,ks) + A{ki + l,k2,k3) 
+A{ki, k2 - 1, fcs) + A{ki, k2 + 1, fca) 

+A{k,, k2, k3-l) + A{ki, k2, ks + 1) (3.28) 

These equations show us the action of H in configurations of the Hilbert space TYg^^ 
for which the three pseudoparticles ( |fci[0]), z = 1,2,3 ) are separated. We already 
know that they are satisfied with the plane wave parametrization ( |3.26| ) and that 

3 

ea = E(-2^i + h + ^7')- (3.29) 

i=i 

For configurations where two pseudoparticles are neighbors at ki and the third pseu- 
doparticles is at A;2 > ^1 + 2, if gives us the following equations 

(es -I- bzi - 22:5 - 2 ^5 -£Zi)A{ki, ki + 1, k2) 
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= A{ki-l,ki + l,k2) + A{ki,ki + 2,k2) 

+A{ki, ki + 1, ^2 - 1) + A{ki, ki + 1, /cs + 1) 

+e fe Bi{ki + 1, fca) + ezeBi{ki, fca) (3.30) 

and a similar set of equations coupling -82(^1,^2) and A{ki, k2, k^), which correspond 
to the meet of two pseudoparticle on the right hand side of the third pseudparticle. 
Comparing the Eq. (|3.30|) with the Eq. (|3.28|) we get a consistency equation 



e zeBi{ki + l,k2) +ezGBi{ki,k2) = A{ki,ki,k2) + A{ki + l,ki + l,k2) 

-Ai A{ki,ki + l,k2) (3.31) 

Similarly, for the right hand side meeting we get 

ezGB2iki,k2 + l)+ezQB2iki,k2) = A{ki, k2, k2) + A{ki, k2 + I, k2 + I) 

-AiA{kuki + l,k2) (3.32) 

These consistency equations are solved by the following parametrization of the 
wavefunctions Bi{ki, k2), i = 1, 2. 

B,{k,,k2) = i?ii(ei6)'^e3' + i?i2(eie3)'^e2^ + 5i3(66)'^et^ 

B2{k,,k2) = 52iei'H6e3)'^+ 522^2^6^3)'^ +^23^3^ (66)'^ (3.33) 

which satisfy the periodic boundary condition provide 

B21 = ^iBi3, B22 = ^2 B12, B23 = Bu (3.34) 
Moreover, the weight functions Bu and -821, i = 1, 2, 3 are determined 

Bu = -Fl2^123 + -^21^213 ; -^21 = -^23^123 + -^32^132 

Bi2 = -^13^132 + -^31^231 5 -^22 = -^13^213 + -^31^312 (3.35) 
-B13 = -^23^312 + -^23^321 ; -^23 = -^12^231 + -^21^321 

where 

Fa, = e^^^^k^-^, a ^6= 1,2,3 (3.36) 

Substituting these relations into the eigenvalue equations (|3.30| ) we obtain the phase 
shift of two pseudoparticles 

^123 ^231 ^ ^132 ^213 ^ ^312 ^123 ^ o^N 

A213 A321 A231 A312 A321 Ai32 
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where 



^ (1 + Ubf - (1 + Ub) {A2^a + As^fc) + ^4Ub + As^g ^/,= 10. 

:i + eae.P-(l + ea6)(A26 + A3ea)+A4ea6 + A5er ^ 

(3.38) 

and the A,, i = 1, 2, 3, 4 are given by ( p.l7|) and (|3.19| ) 



Next, when the three pseudoparticles are neighbors we have the following eigenvalue 
equations 

(es + 4zi - Z5- Z5 -2ezi)A{k, k + l,k + 2) 
= A{k-l,k + l,k + 2) + A{k,k + l,k + 3) 
+e Zq Bi{k + l,k + 2)+ ezeBiik, k + 2) 
+e Zq B2ik, k + 2)+ ezeB2ik, k + 1) (3.39) 

which are automatically satisfied by the above parametrizations. 

In addition to this equations we also have to consider the equations for configu- 
rations where the pseudoparticle of the type \k[—]) and the pseudoparticle \k[0]) are 
separated: 

(eg + Azi - Z5- ^5 - zj -zj)Bi{ki, /ca) 
= B,{k,,k2-l) + B,{k,,k2 + l) 

+Z^Bi{ki - 1, k2) + Z^Bi{ki + 1, k2) 

+ Zq A{ki - 1, ki, k2) + -26^(^1 + 1, fci, k2) (3.40) 

and similar set of eigenvalue equation involving B2{ki,k2), which corresponds to con- 
figurations with the pseudoparticle |A;2[— ]) on the right side hand of the pseudoparticle 
|fci[0]). 

These equations are also satisfied by the above parametrizations. This statement 
was already expected since at this point we always have a far particle as a viewer. 
Therefore, no appeared configurations different from those presented in the sector r = 
2. 

Finally, the action of H on configurations where the two different pseudoparticles 
are neighbors results in two more eigenvalue equations: 

{e3 + 3z,- Z7 -2z^)B,{k,k + l) = Bi{k,k + 2) + B2{k,k + 1) 
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+z^Bi{k - 1, A; + 1)+ A{k - 1, A;, + 1) 

(3.41) 



and 

(e3 + 32l-^7-225)S2(fc,fc + l) = S2(fc-l,fc + l)+5i(fc,fc + l) 

+z:,B2{k, k + 2) + zeAik, k + l,k + 2|3.42) 

Substituting the wavef unctions parametrizations for A{ki, k2, k^) and -Bj(fci, /C2) into 
the equations (|3.41|) and ( |3.42|) and using the relations ( |3.27| ) and ( p.37| ) one can verify 



that they are indeed satisfied. These results tell us that the meeting of the pseu- 
doparticle \k[—]) with the pseudoparticle \k[0]) can be versed as a meeting of three 
pseudoparticle |A;[0]). 

Compounding ( p.37| ) with the periodic boundary conditions ( |3.27|) we arrive to the 



BA equations for the sector r = 3 

3 

^a=I['^ab, a = 1,2,3 (3.43) 

which expresses the factorization of the three pseudoparticle phase shift into the prod- 
uct of two pseudoparticle phase shifts. 

3.5 General sector 

The above results can now be generalized. First we observe that in the sector r > 3 
there is no additional meeting conditions. For example in the sector r = 4 there is a 
meeting of two pseudoparticles of the type |A;[— ]). Nevertheless, we know that the state 
|A;[— ]) is parametrized as two states \k[0]) at the same site and we have verified that the 
meeting of two pseudoparticles ]) behaves as the meeting of four pseudoparticles 
\k[0]) whose phase shift factorizes in a product of two pseudoparticle phase shifts. 

In a generic sector r we build eigenstates of H out of translational invariant prod- 
ucts of Nq one-pseudoparticle eigenstates \k[0]) and two-pseudoparticle eigenstates 
|A;[— ]), such that r = Nq + 2A^_ . These eigenstates are obtained by superposition of 
terms of the form 

!</),) = |0)x|(/>,_i) + |-)x|</),_2) (3.44) 
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with |0o) = 1, |0i) = |0). The corresponding eigenvalue is a sum of single one-particle 
energies 

r 

Er = Nz^ + ^(-2-^1 + h Ha + C') (3.45) 

a=l 

being solutions of the BA equations 

r 

W ^ab, a=l,2,...,r (3.46) 

where 

(i+Ub+ia+^b-{^ + m\ 



ab 



V+Ub + ^a + ib-{^ + '2)ia)' 

a,b ^ 1,2, ...,r. A = 2cosh2?7 (3.47) 

for the ZF model, 

' 1 + - A6 \ A + LCb - {a - ^6 + (A - 2)eaA 



Vl + - A^a A 1 + CaCf " " 6 + (A - 2)^6^ 

a,6 = l,2,...,r. A = 2cosh277 (3.48) 
for the IK model and 

/I + Ub - A6 Wl + Ub + + 6 - (A + 2)^„\ 



\1+Ub-Maj \1+Ub + Ca + Cb-{A + 2)^,J 

a,b = 1,2, ...,r. A = 2cosh2r; (3.49) 



for the osp(l|2) model. 



4 The Algebraic Bethe Ansatz 

In the previous section we have considered the problem of diagonalization of a one- 
dimensional spin chain Hamiltonian using the Coordinate BA. Let us now turn to two 
dimensional classical statistical systems on a lattice. 

Let us consider a regular lattice with N columns and N' rows. A physical state 
on this lattice is defined by the assignment of a state variable to each lattice edge. If 
one takes the horizontal direction as space and the vertical one as time, the transfer 
matrix plays the role of a discrete evolution operator acting on the Hilbert space 
spanned by the row states which are defined by the set of vertical link variables on 
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the same row. Thus, the transfer matrix elements can be understood as the transition 
probabihty of the one row state to project on the consecutive one after a unit of time. 

The main problem now is the diagonalization of the transfer matrix of the lattice 
system. To do this we request the Algebraic BA. 

Again, we start with the graded formulation such that the additional signs are 
represented by e = — 1. Taking e = 1 we recover the non graded cases. 

We recall some basic relations of the graded quantum inverse scattering method. For 
us the basic object will be the R matrix ( |2.6| ), which satisfies -R(0, rj) = pij]) X, where 
Pzf{v) = sinh.rismh.2ri, Pik{v) = — sinhSr] + sinh?7 and Posp(^) = sinh 27] sinh Sr^. 

A quantum integrable system is characterized by monodromy matrix T(A) satisfying 
the equation 



T(A) ® T(/i) 



r(/i) ® T(A) 



(4.1) 



whose consistency is guaranteed by the YB equation (p.5|) . T(A) is a matrix in the 
space V with matrix elements that are operators on the states of the quantum system 
(the quantum space, which will also be the space V). The space V is called auxiliary 
space of T(A). An example of a monodromy matrix is the matrix VR, this follows 
directly from ( |4.1| ). 

The simplest monodromies have become known as C operators, the Lax operator, 
and the monodromy operator T(A) is defined as an ordered product of Lax operators 
on all sites of the lattice: 



T(A) = /:^(A)/:^_i(A)---/:i(A). 



(4.2) 



The Lax operator on the n quantum space is given the graded permutation of ( p.6| 
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/4?(A) 4?(A) 4?(A)\ 



4^(A) (A) LT,>{X) 



(n). 



(4.3) 



V4iHa) 4?(a) 4?(A)y 

Note that Lj-"''(A), i,j = 1, 2, 3 are 3 by 3 matrices acting on the n^^ site of the lattice. 
It means that the monodromy matrix has the form 



T(A) 



Tii(A) ri2(A) ri3(A) 
T2l(A) r22(A) r23(A) 
TsliX) r32(A) r33(A) 



A(A) B,{\) B^{\) 
Ci(A) A2(A) 53(A) 
C2(A) C,{\) A,{\) 



(4.4) 



where 



r.,(A) 



E 4?(A)®i^i::r(A) 

fcl,...,fcjV-l = l 

1,2,3. 



^S-u(A) 



(4.5) 



The vector |0) in the quantum space of the monodromy matrix T(A) that is annihi- 
lated by the operators Tjj(A), i > j (Cfc(A) operators, k = 1,2, 3) and is an eigenvector 
for the operators Tjj(A) ( Ak{X) operators, k = 1,2,3) is called a highest vector of the 
monodromy matrix T(A). 

The transfer matrix r(A) of the corresponding integrable spin model is given by 
the supertrace of the monodromy matrix in the space V, StrT(A). It is the generating 
function of the family of commuting operators in terms of which the Hamiltonian of 
the quantum system is expressed. 

r(A) = StrT(A) = E(-)^(*) T,,(A) = A^{X) + sA^iX) + A,{X) (4.6) 



i=l 



In particular, the Hamiltonians ( |2.12| ) can also be derived by the well-known relation 



^ = «|^(lnr(A)),^o 



(4.7) 



A detailed exposition of the graded quantum inverse scattering method can be 
found in reference p2 . 



In this section we will derive the BA equations of 19-vertex models presented in 



Section 2 using the Algebraic BA developed by Tarasov and recently generalized 
by Martins and Ramos |^ . To do this we need of the commutation relations for entries 
of the monodromy matrix which are derived from the fundamental relation ( |4.1|) . Here 
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these commutation relations do not share a common structure. Therefore, we only 
write some of them in the text and recall ( f4.1| ) to get the remaining ones. 

First of all , let us observe that for each row state one can define the magnon 
number operator which commutes with the transfer matrix of the models 



N 




[r(A),M]=0, M = Y.Mk, = \ 10, (4 

k=l 

This is the analog of the operator used in the previous section and the relation 
between M and the spin total is simply M = N — S^. Once again, the Hilbert 
space can be broken down into sectors In each of these sectors, the transfer 

matrix can be diagonalized independently, r(A)\l/M = m ■ We will now start to 
diagonalize r(A) in every sector: 

4.1 Sector M = 

Let us consider the highest vector of the monodromy matrix T(A) in a lattice of 
sites as the even (bosonic) completely unoccupied state 

^0 = |0) = ®ti j (4.9) 

It is the only state in the sector with M = 0. Using (|4.5|) we can compute the action 
of the matrix elements of T(A) on this reference state: 

A(A)|0) = a;f(A)|0), ^^(A) |0) = a:^(A) |0) , A3(A) |0) = X3^(A) |0) 
Cfc(A)|0) = 0, 5fc(A)|0)^{0,|0)}, A; = 1,2, 3 (4.10) 

Therefore in the sector M = , \I'o is the eigenstate of r(A) = Ai{X) + £^42 (A) + A'i{\) 
with eigenvalue 

Ao(A) = xf (A) + £x^(A) + x^(A) (4.11) 

Here we observe that the action of the operators i?i(A), B2{X) and -B3(A) on the refer- 
ence state will give us new states which lie in sectors with M 7^ 0. 
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4.2 Sector M = 1 

In this sector we have the states Bi{X) |0) and 53(A) |0). Since 53(A) |0) oc Bi{X) |0), 
we seek eigenstate of the form 

^i(Ai)=5i(Ai)|0). (4.12) 

The action of the operator r(A) on this state can be computed with aid of the 
following commutation relations 

A,{X)B^{f,) = z{|,-X)BMA^{X)-^^^^B,{X)A,{^,) (4.13) 
^=<' - '''b.(;.)C3(A) - 4^B,(A)C,(^) (4.15) 



a;3(A — /i) a;3(A — jj,) 

The ratio functions which appear in the commutation relations are defined by 

.(A) = uiX) ' '"^'-^ 



X2(A) ' X3(A)a;4(A) - XQ{X)y^{X) ' 

..(\\ ^3(A) , . X3(A)3;4(A) -X6(A)y6(A) , . 

= RA)- f(-^)=^,(A),.(A)-X3(A)x.(A)- (*-^*" 

When t(A) act on ^i(Ai) , the corresponding eigenvalue equation has two unwanted 
terms: 

r(A)*i(Ai) = (^(A)+£A2(A)+A3(A))^'i(Ai) 

- WA. - A)<(A) ,A) . ||A^.J'(A)1*.(A.) 
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From the matrix elements of 19- vertex models ( p.7| - p.9| ) we can see that X5(A)/x2(A) = 
— ?/5(— A)/x2(— A). Therefore the unwanted terms vanish and \l/i(Ai) is eigenstate of 
r(A) with eigenvalue 



Nfw , 2^('^~'^i) Nf\\ , a;2(A — Ai) ^ 



Ai(A, AO = z{X, - X)xnX) + _ (X) + ^^^^fsW (4-18) 



provided 

/xi(Ai)\ 2 



\X2{Xl 

4.3 Sector M = 2 



= 1 (4.19) 



In the sector M = 2, we encounter two linearly independent states Bi{X)Bi{fi) |0) and 
52(A) |0). (The states B^Bs |0) , B1B3 |0) and B3B1 |0) also lie in the sector M = 2 but 
they are proportional to the state -Bi-Bi |0)). We seek eigenstates in the form 

^2(Ai, A2) = 5i(Ai)fii(A2) |0) + 52(Ai)r(Ai, A2) |0) (4.20) 

where r(Ai, A2) is an operator-valued function which has to be fixed such that ^'2(Ai, A2) 
is unique state in the sector M = 2. 

Here we observe that the operator-valued function r(Ai,A2) is the analog of the 
weight function ^(^1,^2) of the Eq.(p7D. 

It was demonstrated in ||16[ that \l/2(Ai,A2) is unique provided it is ordered in a 
normal way: In general, the operator-valued function \I^„(Ai, A„) is composite of 
normal ordered monomials. A monomial is normally ordered if in it all elements of the 
type Bi{X) are on the left, and all elements of the type Cj{X) on the right of all elements 
of the type Ak{X). Moreover, the elements of one given type having standard ordering: 
Tijj^(Ai)Ti2j2(A2)...Ti„j„(A„). For a given sector M = n, \l/n(Ai, A„) is unique. 

From the commutation relation 

Si(A)Si(/i) = cu(^-A)[5i(/i)5i(A)-— ^— 52(/i)Ai(A)] 

y{fi - A) 

-52(A)Ai(/i) (4.21) 



y(A-/^) 

we can see that ( [4.21 ) will be normally ordered if it satisfies the following swap condition 



^2(A2, Ai) = ^(Ai - A2)^2(Ai, A2) (4.22) 
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This condition fixes r(Ai,A2) in Eg. ( [4.2q ) and the eigenstate of r(A) in the sector 
M = 2 has the form 

*2(Ai, A2) = 5i(Ai)5i(A2) |0) - ^ ■ B2{\i)Ai{\2) |0) . (4.23) 

y[\i - A2) 

The action of transfer matrix on the states of the form ( |4.23| ) is more laborious. In 
addition to ([4.13| - ^A5D and ( [4.21| ) we need appeal to ( [4.1D to derive more eight com- 



mutation relations 

A,{X)BM = ^4^^i?2(/i)A(A)-^4^^i?2(A)Ai(/x) 

-e^^^^MX)B,{f,) (4.24) 
xsin - A) 

A2(A)S2(/i) = z{\-fl)z{fI-\)B2{fi)A2{X) 



a;2(A-/i) X2(A-/i)' 
A.{\)BM = ^4^52(/x)A3(A) - M^5,(A)A3(/.) 

XsiX-IJ,) X3(A-/i) 



(4.25) 



y(A - /i) 



^ B3{X)BM (4.26) 



Ci(A)i?i(/i) = eB,{f,)C,{X) + ^\ [Ai(/i)A2(A) - A,{X)AM] (4.27) 

X2(A - fij 

C3(A)Si(/i) = ^ '''!^~^| i?iMC3(A)- '^^|^~^| i?i(A)C3(/i) 

X3(A-/i) X3(A-/i) 

■ ^ -[Ai(^)A3(A) - A2{X)A2{fl)] + ^^^B2{fl)C2{X) 



y{X-n) x^iX-jj)' 

(4.28) 

BiiX)BM = -^^^f2{l^)Bi{X) + ^'|^I^) ^iM^2(A) (4.29) 

5i(A)i?3(/i) = ei?3(/i)i?i(A)-^4T^^2(/x)A2(A) + ^4T^^2(A)A2(^) 

X2(A-/i) a;2(A-yu) 

(4.30) 

52(A)Si(^) = -^^^fMB2{X) + ^^1^ ~ |^| i?2(^)i?i(A) (4.31) 

Here we observe that in this approach the final action of r(A) on normally ordered states 
must be normal ordered. This implies in an increasing use of commutation relations 
needed for the diagonalization of t(A). For example, the action of the operator Ai{X) 
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on ^'2(Ai, A2) has the form 



A(A)*2(Ai,A2) = ^(Aio)^(A2o)xf(A) *2(Ai,A2) 

-^(A2i)xf (Ai) Si(A)Si(A2) |0) 



X5(A2o) -(Ai2U^^^^ B,{X)B,{X,)\0) 



a^2(A2o) a;(Ai2) 
_^ / ^(Aio) X5(A2o) 1 _^ a:;7(Aio) 



^(Aio) a:^2(A2o) y(Aoi) a;3(Aio) y(Ai2) / 
xxf (Ai)xf (A2) B^iX) |0) (4.32) 

where Aab = Xa — Xb, a ^ b = 0, 1,2, with Ao = A. Here we have used the following 
identities satisfied by the matrix elements of these 19-vertex models: 

z{Xab) X5{Xcb) ^ XQ^Xab) 1 _ X^jXqb) X^^Xca) _^ z{Xac) X^^Xcb) 



Uj{Xab) X2{Xcb) XsiXab) y(Xac) X2{Xab) X2{Xca) <^(Aac) 2^2 (Acft) 

.x^iXcb) 1 , xi{Xab) 1 r\ \ r\ \ 

z{Xab) /X ^ M ^ + f\ \ f\ \ ^ 4Kb)z{Xcb^ 



X2{Xcb) y{Xab) XsiXab) yiXac) ViKc) 

Uj{XabMXba) = 1, (qt^Mc) (4.33) 

Similarly, for the operator A2(A) we have 

^2(A)*2(Ai, A2) = £^4r44r4^^ W 

-£^^4v4z(A2i)x^(A2) i?i(A)i?i(AO |0) 



_2y5(Aoi) 2:(Ai2) N 



X2(Aoi) a;(Ai2; 



x^{Xi) 5i(A)5i(A2)|0) 



+z{^2i)^^x^iXi) Bs{X)B,{X2) |0) 

i?3(A)Si(A0 |0) 

'^(Al2) |/(Ao2) 

^ ^2 y5(Aoi) / y5(A2l) 1 _^ ^(Aqi) 1 _ y5(Aoi) 1 \ 
a^2(Aoi) V^2(A2i) y(Xoi) u;{Xoi) y(Ao2) a;2(Aoi) y(Ai2) / 
xxf (A2)x^(Ai) S2(A) |0) 

In this case we have used more two identities: 

z{Xab) 1 ^ 1/5 (Acb) 1 _ ?/5(Aab) 1 ^ 2:(Afec) 1 



•^(Aafe) y(Aac) a;2(Ac6) y(A„;,) a:;2(Aa6) y(A6c) "^(Aftc) y(A„c) 
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cb I 



ac) ^ VbiKb) VbiKc) 



X2{Kc) X2{Xab) X2{Xbc) X2{Xa 

a 7^ b c 



(4.35) 



Finally, for A3 (A) we get 



A3(A)^2(Al,A2) 



a^2(Aoi) X2(Ao2) 

3^3(^01) 3;3(Ao2, 



x^(A) ^2(Ai,A2 



12) 



+ 



t^(Ai2) |/(Aoi) 
y7(Aoi) 1 2/5(Aoi 



X^^(Ai) i?3(A)i?i(A2) |0)-£^z(A2l) 



y(A 



<(A2) B3(A)5i(Ai)|0) 



02 J 



4nAi)x^(A2) 52(A) |0) 



(4.36) 



xsiXoi) yiXu) xs{Xoi)y{Xo2) 

Here we also have used the identities ( |4.33| ) and ( [4.35| ). 

From these relations one can see that all unwanted terms of r(A)^2(Ai, A2) vanish. 
It means that \E'2(Ai, A2) is an eigenstate of the transfer matrix r(A) with eigenvalue 



,N/\\ , ,3 ^(Aoi) ^(Ao2) N/w I 3;2(Aoi) a;2(Ao2) at. 



A2(A, Ai, A2) = z(Aio)^(A2o)xf (A) + ^^^^^^(A) + <(A) 

uj{Xoi) uj{Xo2) a;3(Aoi) a;3(Ao2) 

(4.37) 

provided the rapidities Ai and A2 satisfy the BA equations 

^£liM]".,3£iM^A.,. (4.38) 

Va;2(Aa)/ z{Xba) 

4.4 General Sector 

The generalization of the above results to sectors with more than two particles proceeds 
through the factorization properties of the higher order phase shifts discussed in the 
previous section. Therefore, at this point we shall present the general result: In a 
generic sector M = n , we have n — 1 swap conditions 

^n(Ai, ■ ■ ■ , Aj-i, Aj+i, Aj, ■ ■ ■ , A„) = uj{Xi — Ai+i)\E'n(Ai, ■ ■ ■ , Ai_i, A,, Aj+i, ■ ■ ■ , A„) 

(4.39) 

which yield the n — 1 operator-valued functions rj(Ai, ■ ■ ■ , A„) . The corresponding 
normal ordered state \E'„(Ai, ■ ■ ■ , A„) can be written with aid of a recurrence formula 
1§: 

^„(Ai,...,A„) = $„(Ai,...,A„)|0) (4.40) 
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where 



$„(Ai, A„) = Si(Ai)$„_i(A2, A„) 

n 1 n ^ 

-^2(Ai)E-7T TT n 2(A,-A,)$„_2(A2,...,A„-,A„Mi(A,) (4.41) 

j=2 yy^i - ^j) k=2,k^3 

with the initial condition — $i(A) = 5i(A). 
The scalar function Z{\k — \j) is defined by 

^(A. - A,) = I 1^ " ^1 I l"^'- (4.42) 

The action of the operators Aj(A), i = 1, 2, 3 on the operators have the following 
normal ordered form 



Ai(A)$„(Ai, A„) = n •^(Afc - A)$n(Ai, A„)Ai(A) 

fc=i 

-^iWE?#^ n ^(A.-A,)$„_i(Ai,...,A„...,A„Mi(A,) 

n j — 1 n 

+^2(A) J2 E CiKA, A,, A,) n ^(Afe - Ai)2^(Afe - A,) 

j=2 Z=l k=l,k^j,l 

x$„_2(Ai, Ai, Aj-, Xn)Ai{Xi)Ai{Xj) (4.43) 
where Gji{X, A;, Aj) are scalar functions defined by 

r (\ \ \ 1 = ~ ^ I ~ ^^^^^ ~ ^ a 44^ 

Uji^A, Au A,) ^^^^^ _ ^^^^ _ + ^^^^ _ ^^^^^ _ _ ^4.44J 

For the action of A^{X) we have a similar expression 

yl3(A)$n(Ai, AO = n ^4^^^n(Ai, A0A3(A) 

fe=l ^Sl'^ — Ak) 

-£-S3(A)E— — — n 2(A,-A,)<fn-i(Ai,...,A„...,A0A2(A,) 

j=l k=l,k^j 

n j — 1 n 

+^2(A) ^ Yl Hji{X, Xi, A,) n ^(A.- - Afe)Z(A; - A^) 

j=2 1=1 k=l,k^j,l 

x$„_2(Ai, Xi, Ai, A„)A2(A0A2(A,) (4.45) 
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where the scalar functions Hji{X, Xi, Xj) are given by 

^" - X3(A-A0y(A.-A,) ' xsiX-XOviX-Xj) ^^"^^^ 
The action of the operator A2{X) is more cumbersome 

A2(A)<l>„(Ai, A„) = 5" n \. X > n(Ai, A0A2(A) 

fe=i ^(A - Afc) 

-£'^Si(A)X:^4t-TT ft ^(A,-A,)$„_i(Ai,...,A„...,A0A2(A,) 

" 1 " A 

+^3(A)E a X ^ n ^(Afc-A,)«l>„_i(Ai,...,A„-,An)Ai(A,) 

{n j—l n 
^ ^ F,,(A, A;, A,) n ^(-^fc - -^O^lAi - Afc) X 
j=2 1=1 k=l,k^j,l 

$„_2(Ai, A;, Aj, Xn)Ai{Xi)A2{Xj) + 

Ti i — 1 n 

5] 5] F,,(A, A,, A,) n ^(Ai - Afe)Z(Afe - A,) x 

3=2 1=1 k=l,k=^j,l 

*„_2(Ai, Xi, A,-, A„)Ai(A,-)^2(A0| (4.47) 
where we have more two scalar functions 



F(XX \ ) = 1/5(A-A;) j yrMi-^j) 1 I ^(A-Af) 1 
''^ ' a;2(A-A0 lx2(A,-A,)y(A-A0 ^A-AOz/(A-A,) 



a;2(A - A;) y{Xi - Xj) 



(4.48) 



Y (\ \ A ^ - ^ f A^ ^5(A-A,) ?/5(A-A0 i/5(A;-A,) | 

- ^(A^ 1^^^ - ^^^^5a^ " X2(A-A0x2(A,-A,)l ^'-''^ 
Prom these relations immediately follows that ^m{Xi, Am) are the eigenstates of 
r(A) with eigenvalues 

Am = x,(A)^ n .(A„-A)+s^+^a;2(A)^ f] 4^^+^3(A)^ fl "f.'^A (4-50) 

a=l a=l - ^a) a=l ^^[^ " '^aj 

provided their rapidities A^, i = 1, M satisfy the BA equations 

fxi^\\^M-,i fi^^^(A,-A„), a = l,2,...,M (4.51) 
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To conclude this section we remark that equations ( [4.50| ) and ( [4.51| ) reproduce the 
known results in the literature: Using ( |2.(j| ) we reproduce the BA solution of the IK 
model |T6[ and using ( p.7| ) they are the BA solution for the ZF model derived by a fusion 



procedure in [|13|, [|14| and . Specifically, in the case of the rational solution and for 
e = —1 we obtain previous results derived by the Analytical and Algebraic ||2^ BA 
approach for the rational osp{l\2) vertex model. Furthermore, by using the equation 
( [4 .7]) we recover the expressions derived in the previous section via the Coordinate BA. 



5 Conclusion 



In the first part of this paper we have applied the Coordinate BA to find the spectra 
of Hamiltonians associated to three 19-vertex models, including a graded model. This 
procedure was carried out for periodic boundary conditions. 

We believe that the method here presented could also be applied for Hamiltonians 
associated with higher states vertex-models. For instance, in the quantum spin chain 



s = 3/2 XXZ model we have four states: 



and 



. It means 



that the state |A;(l/2)) can be parametrized by plane wave and the states and 



kW\) at the same site , respectively, multiplied by some 



^[ 2 V ^^"-^ three states 

weight functions. 

These weight functions are responsible by the factorized form of the phase shift 
of two particle ( p.38| ). In the ZF model we do not a factored form for the two- 
pseudoparticle phase shift because its weight function ( p.l7| ) is a constant. It means 
that the state ]) behaves exactly as two states \k\Q\) at the same site. This is in 
agreement with the fact that the ZF model can be constructed by a fusion procedure 
of two six-vertex model. 

In the second part of this paper we have applied the Algebraic BA to find the 
spectra of the transfer matrices of these three state vertex models. The method here 



presented was developed by Tarasov [|T^] and generalized by Martins ||21|, |2^. It is 
general enough to include the ZF model as well as the graded osp{l\2) model. 

There are several issues left for future works. A natural extension of this work is to 
consider these Bethe Ansatze with open boundary conditions via refiection matrices. 
The transfer matrix of the ZF model with the most general diagonal refiection matrix 
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has been diagonalized by Mezincescu at al by generalizing the fusion approach used 
to solve the corresponding model with periodic boundaries. BA equations for both the 
ZF model and IK model with open boundaries were derived by Yung and Batchelor in 



ref.[^. Nevertheless, basead on the Tarasov-Martins approach, the Algebraic BA of 
the IK model with a diagonal i^'-matrix was recently derived by Fan ||26|| . 
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LIST OF CHANGES 



Reference: A/99391/PAP 

Title: Bethe Ansatze for nineteen- vertex models 
Author: A. Lima-Santos 

1. page 27, before the section 5 (Conclusion) we replace all text (hues 1,2 and 3): 
In particular, via the Coordinate BA. by the following text: 

To conclude this section we renicirk that equations (4.50) and (4.51) 
reproduce the known results in the literature: Using (2.8) we reproduce 
the BA solution of the IK model [16] and using (2.7) they are the BA 
solution for the ZF model derived by a fusion procedure in [13], [14] and 
[15] . Specifically, in the case of the rational solution and for s = —1 we 
obtain previous results derived by the Analytical [20] and Algebraic [23] 
BA approach for the rational osp(l|2) vertex model. Furthermore, by using 
the equation (4.7) we recover the expressions derived in the previous section 
via the Coordinate BA. 
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